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Abstract 

We performed a combined analysis of the parameter space of the Minimal 
Supersymmetric Standard Model (MSSM) taking into account cosmological 
and accelerator constraints including those from the radiative b — > 57 decay 
measured by the CLEO collaboration. 

Special attention is paid to the event rate, R, of direct dark matter 
neutralino detection. We have found domains of the parameter space with 
R ~ 5 — 10 events/kg/day. This would be within the reach of current dark 
matter experiments. The b — > 37 data do not essentially reduce these large 
event rate domains of the MSSM parameter space. 
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1 Introduction 



Rare b — > s'j decay observed by the CLEO collaboration |L| with the branch- 
ing ratio Br(b —>■ sj) = (2.32±0.67) ■ 10 -4 , has been recognized as a stringent 
restriction for physics beyond the standard model (SM). 

In the Minimal Supersymmetric Standard Model (MSSM) || this decay 
proceeds through 1-loop diagrams involving W-boson, charged Higgs boson, 
chargino, neutralino and gluino || |J. Since the SM prediction is consistent 
with the measured branching ratio BR(b —>■ 37) the MSSM contributions are 
stringently restricted. A dramatic reduction of the allowed MSSM parameter 
space due to the b — > 37 constraint was reported by many authors [|5], The 
impact of this constraint on prospects for direct detection of the dark matter 
(DM) neutralino (x) via elastic scattering off various nuclei has been also 
analyzed. In || []] it was found that within popular supergravity models the 
detection rate becomes too small for observation if the CLEO constraint is 
incorporated in the analysis. In this case the upper bound on BR(b — > 57) 
implies a stringent lower bound on the mass of the pseudoscalar Higgs boson 
(7714.) of the MSSM if sparticles are heavy. It leads to strong suppression of 
the elastic neutralino-nucleus scattering cross section. 

Scenarios with lighter sparticles and pseudoscalar Higgs have been studied 
in the literature as well ||. Also here, first results for direct detection of 
neutralinos were pessimistic. Again, only small counting rates had been 
found in the domain of the MSSM parameter space satisfying the BR(b — > 
57) constraint. 

A more comprehensive exploration of the constrained MSSM parameter 
space f7|, ||, U discovered, however, that the b — > 57 constraint makes actually 
only a moderate effect on the expected event rate. It was realized that 
relatively light neutral Higgs bosons, leading to typically large event rates, 
are compatible with this constraint in contrast with the results of previous 
analyses. A key observation is that the charged Higgs boson contribution to 
BR(b — > 57) can be efficiently compensated by the chargino contribution in 
a large domain of the MSSM parameter space. This essentially relaxes the 
lower bound both for charged and for neutral Higgses imposed by the CLEO 
result. 

Recently an exciting result was obtained in A sophisticated scan of 



the MSSM parameter space constrained by the known experimental bounds 
including BR{b — > s-y) picked up points with an unexpectedly large detection 
rate of DM neutralinos with a mass around 1 TeV. For a germanium ( 76 Ge) 
target an integrated detection rate was found at a level of 10 events/kg/day 
and even up to 100 events/kg/day for sodium iodide (Nal). 
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If these results are correct they have important consequences for direct 
DM detection. Such large allowed event rates would mean that the cur- 
rent DM experiments have already entered an unexplored part of the MSSM 
parameter space. 

Certainly, before such a conclusion can be made one has to be sure that 
the above cited results, obtained in a specially arranged scan, are not an 
artifact having no relation to physics. One may suspect, for instance, a 
specific instability of the numerical code used in the analysis. Note, that 
the standard scan without special sampling of the model parameters did not 
produce even a single point with an event rate larger than 1 event/kg/day. 
Therefore, an independent search for large event rate points within the MSSM 
parameter space is apparently demanded. 

In the present paper we carry out a systematic scan of the MSSM param- 
eter space constrained by the known accelerator data and by the requirement 
that the DM neutralinos do not overdose the universe. We adopt the unifica- 
tion scenario [11 with a non-universal scalar mass when the soft Higgs mass 
parameters are not equal to the common sfermion soft mass parameter at the 
unification scale. In this case the Higgs and sfermion masses are not strongly 



correlated parameters. As discussed in [ XT | , this minimal relaxation of the 
complete unification conditions allows one to avoid one of the most stringent 
theoretical limitations on the allowed values of the neutralino detection event 
rate. Other unification conditions do not make such an effect and tolerate 
large event rate values. Therefore, we keep these unification conditions to 
reduce a number of free parameters. The latter is crucial for a fine scanning 
of the MSSM parameter space which we are going to carry out. 

The aim of the scan is to "detect" those domains in the parameter space 
where the event rate R of the direct DM detection approaches experimentally 
interesting values R >1 events/kg/day. Applying an extensive standard scan 
procedure we have found such domains with a detection rate of about 10 
events/kg/day in 73 Ge. We demonstrate that incorporation of the b — > 37 
constraint leads to only a moderate effect on these domains. 



The special sampling described in |10[ has also been applied. We did 
not reproduce large event rate domains located around a value of 1 TeV 
for the neutralino mass, reported in the cited paper. No points with R > 1 
event /kg/ day have been " detected" in this region of neutralino masses neither 
in the standard nor in the special scans. 

The paper is organized as follows. In section 2 we specify the MSSM and 
give a list of the formulas relevant to our analysis. In section 3 we summarize 
the experimental inputs for our analysis and in section 4 we summarize the 
formulas for event rate calculations, then in section 5 discuss our numerical 
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procedure and results. Section 6 contains the conclusion. 



2 Minimal Supersymmetric Standard Model 

The MSSM is completely specified by the standard SU(3)x SU(2)x U(l) 
gauge couplings as well as by the low-energy superpotential and "soft" SUSY 
breaking terms ||. The most general gauge invariant form of the R-parity 
conserving superpotential is 

W = h E UE c E\t l3 + h D Q j D c Hie i:j + huQ j U c H l 2 e i:j + ////;//^, ; (1) 

(ei2 = +1)- The following notations are used for the quark Q (3,2,1/6), 
D c (3,1,1/3), U c (3,1,-2/3), lepton L (1,2,-1/2), E c (1,1,1) and Higgs 
E x (1,2,-1/2), H 2 (1,2,1/2) chiral superfields with the SU(3) c xSU(2) L x 
U(l)y assignment given in brackets. Yukawa coupling constants h E) D,u are 
matrices in the generation space, non-diagonal in the general case. For sim- 
plicity we suppressed generation indices. 

In general, the "soft" SUSY breaking terms are given by [I2|: 

C SB = -WM A \A\A-m^H^-ml 2 \H 2 \ 2 -ml\Q\ 2 

1 A 

- m%\Df - m 2 \U c \ 2 - m\\L\ 2 - m%\Ef 

- (h E A E VE c Hle tJ + h D A D QW c Hie tJ 

+ h u A u Q'U c H t 2 e ij + h.c) - i /)'/'//; //2'u + h.c) (2) 

As usual, M 3j2 ,i are the masses of the SU (3) x SU (2) x U(l) gauginos g, W, B 
and rrii are the masses of scalar fields. Al, Ad, Ay and B are trilinear and 
bilinear couplings. 

Observable quantities can be calculated in terms of the gauge and the 
Yukawa coupling constants as well as the soft SUSY breaking parameters and 
the Higgs mass parameter [i introduced in Eqs. Under the renormal- 

ization they depend on the energy scale Q according to the renormalization 
group equations (RGE). 

It is a common practice to implement the grand unification (GUT) con- 
ditions at the GUT scale Mx- It allows one to reduce the number of free pa- 
rameters of the MSSM. As explained in the introduction, we adopt a scenario 
with a non-universal Higgs mass with the following set of GUT conditions: 

m L (M x ) = m E (M x ) = m^Mx) = m v (M x ) = m D (M x ) = m , 
m Hl (M x ) = m H2 (M x ), (3) 
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Au(M x ) = A D (M X ) = A L (M X ) = A , 

Mi{M x ) = m V2 , (4) 

5 g' 2 g 2 g 2 

ai(M x ) = a GUT , where a x = -—, a 2 = — , a 3 = -p-, (5) 

3 47r Air An 

g' , g and g s are the U(l), SU(2) and SU(3) gauge coupling constants. As seen 
from Eqs.([|) and (|3]) the Higss mass parameters mj/ 12 are not equal to the 
common sfermion mass mo at the GUT scale M x . Accepting the GUT con- 
ditions above, we end up with the following free MSSM parameters: the com- 
mon gauge coupling acuT', the matrices of the Yukawa couplings hf 1 , where 
i — E,U, D\ soft supersymmetry breaking parameters mo, m 1 / 2 , A Q , B, the 
Higgs field mixing parameter fi and an additional parameter of the Higgs 
sector uia being the mass of the CP-odd neutral Higgs boson. Since the 
masses of the third generation are much larger than masses of the first two 
ones, we consider only the Yukawa coupling of the third generation and drop 
the indices a, b. 

Additional constraints follow from the minimization conditions of the 
scalar Higgs potential. Using these conditions the bilinear coupling B can 
be replaced in the given list of free parameters by the ratio tan/? = Vz/vi of 
the vacuum expectation values of the two Higgs doublets. 

We calculate the Fermi-scale parameters in Eqs.(|l|) and (0) in terms of 
the above listed free parameters on the basis of 2-loop RGEs following the 
iteration algorithm described in |TB . 



The Higgs potential V including the one-loop corrections AV can be 
written as: 



V(H°,H°) = ml\H°\ 2 + ml\H°\ 2 -ml(H° 1 H° + h.c.) 

2 '2 

+ ^4^-(|tfiT-|tf 2 T) 2 + A^, (6) 



with AV = ^^(-1) 2 ^(2J, + 1)C^ 



' m 2 3 
V"2 



where the sum is taken over all possible particles with the spin Jj and with 
the color degrees of freedom Cj. The mass parameters of the potential are 
introduced in the usual way as 

m l,2 = m H li2 + /A m l = ( 7 ) 

They are running parameters with the scale Q-dependence m^Q) determined 
by the RGE. The 1-loop potential (0) itself is Q-independent up to, field- 
independent term depending on Q, irrelevant for the symmetry breaking. 
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At the minimum of this potential the neutral components of the Higgs 
field acquire non-zero vacuum expectation values (H® 2 ) = i> 1)2 triggering the 
electroweak symmetry breaking with g 2 (vf + u|) 

The minimization conditions read 



2m\ 
2ml 



2ml tan f3~ M z cos 2/5 - 2S i 
2ml cot /3 + Ml cos 2/3 - 2£ 2 , 



(8) 
(9) 



where 



fr> with ^ 



HeHf 2 , are the one- loop corrections fL4| : 



32tt- 



E(-1) 2Ji (2^ + 1) 



1 (9m? 




(10) 



As remnant of two Higgs doublets ifi )2 after the electroweak symmetry break- 
ing there occur five physical Higgs particles: CP-odd neutral Higgs boson A, 
CP-even neutral Higgs bosons H, h and a pair of charged Higgses H ± . Their 
masses m A , m h H , m H ± can be calculated including all 1-loop corrections 
as second derivatives of the Higgs potential in Eq. (0) with respect to the 
corresponding fields evaluated at the minimum jTS], |TBJ. 

The neutralino mass matrix written in the basis (B, W 3 , H®, H®) has the 



form: 



M, 



( M l 


-M z cps % 
V M z sps w 





M 2 
M z c p c w 
-M z s p c v 



-M Z C(3S h 

M z c p c w 




M z sps w 

-M z Sf3C v 

-H 




;n) 



where s w = sin^, c w = cosOw and Sp = sin/5, eg = cos/5. Diagonalizing 
the mass matrix above by virtue of the orthogonal matrix M one can obtain 
the four physical neutralinos x% with the field content 



Xi 



Aii-B + N i2 W 3 + JVai?? + M iA Hl 



(12) 



and with masses m Xi being eigenvalues of the mass matrix (|TT| ) . The lightest 
neutralino Xi we denote x- I n our analysis x is the lightest SUSY particle 
(LSP). 

The chargino mass term is 




+ h.c. 



with the mass matrix 



M 2 V2M W sin (3 

\/2M w cos (3 fj, 



(13) 



(14) 
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which can be diagonalized by the transformation 

X- = U il W- + U i2 H~, x + = V a W + + V t2 H+ (15) 
with U*M.^±V^ = diag(M~±, M~±), where the chargino masses are 

ML =l\Ml + fi 2 + 2M* v T 



{Ml - fi 2 ) 2 + 4M^ cos 2 2(3 + 4M^(M| + fi 2 + 2M 2 /i sin 2(3) 

The mass matrices for the 3-generation sfermions t, b and f in the /l-/r 
basis are: 



.A/f 2 = 

* I m t (A t - [i cot /?) m 2 7 + m 2 - § (M^ - M|) cos 2/3 



mk+m 2 t + \{AM$r - M§) cos 2/? m t (A t - /xcot /?) 



Ml = 



m 2 Q + m 2 b - \{2M^ + M 2 Z ) cos 2(3 m b {A b - fj, tan (3) 

m b (A b - /itan/3) m 2 ^ + mf + ± (M^ - Mf) cos 2(3 

m\ + m 2 - \{2M^ - Mf)cos2/3 m T (A T - fi tan (3) \ 

m T (A T - jj, tan (3) m| + m 2 + (M^ - Mf ) cos 2/3 J 



For simplicity we ignored in the sfermion mass matrices a non-diagonality 
in the generation space which is important only for the b — > S7 -decay. 

3 Constrained MSSM parameter space 

In this section we shortly summarize the theoretical and experimental con- 
straints used in our analysis. 

Solution of the gauge coupling constants unification (see Eq.(D)) using 
2-loop RGEs allows us to define the unification scale Mx- The following 
standard definitions are used: a,\ = 5a/ (3 cos 2 W ), a 2 = a/ sin 2 6^. The 
world averaged values of the gauge couplings at the Z° energy were obtained 
from a fit to the LEP data |I7), M w JTJ and m t [Lj, g§: oT x (M z ) = 
128.0 ± 0.1, sin 2 ^ = 0.2319 ± 0.0004, a 3 = 0.125 ± 0.005. The value of 
a~ 1 (Mz) was updated from |2TJ by using new data on the hadronic vacuum 
polarization ||22|| . 

SUSY particles have not been found so far and from the searches at LEP 
one knows that the lower limit on the charged sleptons is half the Z° mass 
(45 GeV) to be above 60 GeV For the charginos the preliminary lower limit of 
65 GeV from the LEP 140 GeV run was used The lower limit on the lightest 
neutralino is 18.4 GeV while the sneutrinos have to be above 41 GeV 
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Radiative corrections trigger spontaneous symmetry breaking in the elec- 
troweak sector. In this case the Higgs potential has its minimum for non-zero 
vacuum expectation values of the fields. Solving for Mz from Eqs. (|8|) and 
©yields: 

Mz 2 _ m\ + T, 1 - (mj + S 2 ) tan 2 (3 
2 tan 2 /3-l ' 1 J 

where the Hi and E 2 are defined in Eq. (|T0|) . This is an important constraint 
which relates the true vacuum to the physical Z-boson mass M z = 91.187 ± 
0.007GeV. 

Another stringent constraint is imposed by the branching ratio BR(b — > 
57) measured by the CLEO collaboration [Q to be: BR{b — > sj) = (2.32 ± 
0.67) x 10" 4 . 

In the MSSM this flavour changing neutral current (FCNC) receives in 
addition to the SM W — t loop contributions from H ± — t, x — t and g — q 
loops. The x ~~ t loops, which are expected to be much smaller, have been 
neglected [0, |26]| . The g — q loops are proportional to tan/3. It was found 
[[0| that this contribution should be small, even in the case of large tan/? 



and therefore it was neglected. The chargino contribution, which becomes 
large for large tan/? and small chargino masses, depends sensitively on the 
splitting of the two stop masses. 

Within the MSSM the following ratio has been calculated M: 



BR{b -> s 7 ) 
BR{b -»• cev) ~ 

I V t *Vtb 1 2 r, Q c D Ga 



(17) 



rj^A + f (r/ 14 / 23 - V 16 / 23 )A g + C 



\V cb \ 2 L it I(m c /m b )[l - (2/3n)a s (m b )f(m c /m b )Y 
where 

C w 0.175, / = 0.4847, 77 = a s (M w )/a s (m b ), f(m c /m b ) = 2.41. 

Here f(m c /m b ) represents corrections from leading order QCD to the known 
semileptonic b — > cev decay rate, while the ratio of masses of c- and b- 
quarks is taken to be m c /m b = 0.316. The ratio of CKM matrix elements 
^y V ^2 — 0.95 was taken from ||27|| , the next leading order QCD-corrections 
from |28[ . A„ „ are the coefficients of the effective operators for fe-7 and for 



bs-g interactions respectively. 

Assuming that the neutralinos form a dominant part of the DM in the 
universe one obtains a cosmological constraint on the neutralino relic density. 
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The present lifetime of the universe is at least 10 10 years, which implies 
an upper limit on the expansion rate and correspondingly on the total relic 
abundance. Assuming h > 0.4 one finds that the contribution of each relic 
particle species x has to obey 



n x h 2 < 1, 



(18) 



where the relic density parameter Q x = p x j p c is the ratio of the relic neu- 
tralino mass density p x to the critical one p c = 1.88 ■ 10~ 29 liQg-cm~ 3 . 

We calculate fl x hl following the standard procedure on the basis of the 



approximate formula [31, 32 



fl x hl 



2.13 /T v 



10 



ii 



2.7 K° 



N 



1/2 



GeV- 



axp + bxp/2 



(19) 



Here T 7 is the present day photon temperature, T x /T 7 is the reheating factor, 
xp = T F /m x « 1/20, T F is the neutralino freeze-out temperature, and N F 
is the total number of degrees of freedom at Tp. The coefficients a, b are 
determined from the non-relativistic expansion 



< o- annyh v a + bx 



(20) 



of the thermally averaged cross section of neutralino annihilation. We adopt 
an approximate treatment not taking into account complications, which oc- 
cur when the expansion ( |20| ) fails j3(|. We take into account all possible 
channels of the x ~ X annihilation. The most complete list of the relevant 
formulas for the coefficients a, b and numerical values for the other parameters 

can be found in |32 . 



in eqs. ([LJ) and (g 

Since the neutralinos are mixtures of gauginos and higgsinos, the annihi- 
lation can occur both, via s-channel exchange of the Z° and Higgs bosons and 
t-channel exchange of a scalar particle, like a selectron ]3"3|]. This constrains 



the parameter space, as discussed by many groups |7|, |34], |35|. The size 
of the Higgsino component depends on the relative sizes of the elements in 
the mixing matrix (PD, especially on tan/? and the size of the parameter p. 

In the analysis we ignore possible rescaling of the local neutralino density 
p which may occur in the region of the MSSM parameter space where Q x h 2 < 



0.025 P, [3(| |37|. This is a minimal value corresponding to DM concentrated 
in galactic halos averaged over the universe. If the neutralino is accepted as a 
dominant part of the DM its density has to exceed the quoted limiting value 
0.025. Otherwise the presence of additional DM components should be taken 
into account, for instance, by the mentioned rescaling ansatz. However, the 
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halo density is known to be very uncertain. Its actual value can be one order 
of magnitude smaller. Therefore, one can expect that the rescaling takes 
place in a small domain of the MSSM parameter space. Another point is 
that the SUSY solution of the DM problem with such low neutralino density 
becomes questionable. We assume neutralinos to be a dominant component 
of the DM halo of our galaxy with a density p x = 0.3 GeV-cm -3 in the solar 
vicinity and disregard in the analysis points with fl x h 2 < 0.025. 



4 Neutralino- Nucleus Elastic Scattering 

A dark matter event is elastic scattering of a DM neutralino from a target 
nucleus producing a nuclear recoil which can be detected by a suitable de- 
tector. The corresponding event rate depends on the distribution of the DM 
neutralinos in the solar vicinity and the cross section a e i(xA) of neutralino- 
nucleus elastic scattering. In order to calculate a e i(xA) one should specify 
neutralino-quark interactions. The relevant low-energy effective Lagrangian 
can be written in a general form as 

L eff = E (A • WKX " qYlrt + j± ■ C q ■ XX ■ qq) + O (j^j , (21) 

where terms with vector and pseudoscalar quark currents are omitted being 
negligible in the case of non-relativistic DM neutralinos with typical velocities 
v x pa 10~ 3 c. 

In the Lagrangian ( |2"T| ) we also neglect terms which appear in supersym- 
metric models at the order of l/m~ and higher, where m q is the mass of the 
scalar superpartner q of the quark q. These terms, as recently pointed out 
in P2 |, are potentially important in the spin- independent neutralino-nucleon 
scattering, especially in domains of the MSSM parameter space where m q is 
close to the neutralino mass m x . Below we adopt the approximate treatment 



of these terms proposed in [32] which allows "effectively" absorbing them 



into the coefficients C q in a wide region of the SUSY model parameter space. 

nl r \f?. - A/?. 
M/(V L 2 

cos 2 6 q 4> 2 qL + sin 2 e q 4> 2 qR ) 

2/3 j.2 , 2 a j.2 





14 ' JLO rp 




L 2 








(m x + m q ) 2 







m h - ( m x + m i) 



sin 9 q (f) qL + cos 9 q (j) qR ) 
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+ 



(m x + m q ) 2 m~ 2 — (m x + m q y 



M w P q sin 25, T 3 (AT 12 - ttme w N n ) 



m\\ - (m x + m q ) 2 m 2 q2 - (m x + m q f 



(22) 



C a 



Here 



2 1 
m h 

^cos 2 



77 1 



H 



sin 6L 



cos 2 0„ 



- sin 2 9. 



m\ x - [m x 



+ sin 26 ' 



4Mvk 
1 



-Pi 



+ m q ) 

Mw L 



m 2 q2 ~ 



(m x + m q ) 



q YqL 
2 



>qL (PqR) 



{m x + m q ) 



m x + m q ) 



F H 

h q 

Hq 

<\>qL 
4>qR 



(A/12 — A/ii tan^)(A/i4 cos«h +A/13 sina^), 
(A/12 — A/ii tan 0^) (A/14 sin a# — A/13 cos , 



. 1 , COS «;? 

2 +T3) ^" 



1 sin 

2 cosp 



. 1 sin an , 1 „ cos a# 
2 smp 2 cosp 



tan^ Q A/11, 
2 smp 



T<i) tan 8 



w ■ 



13 



cos /3 



(23) 



Our formulas for the coefficients A q and C q of the effective Lagrangian take 
into account squark mixing q~L — q~R and the contribution of both CP-even 



Higgs bosons h, H. The formulas coincide with the relevant formulas in |p2 
neglecting the terms ~ l/ 777 ! and higher. These terms are taken into account 
"effectively" by introducing an "effective" stop quark t propagator. 

A general representation of the differential cross section of neutralino- 
nucleus scattering can be given in terms of three spin-dependent Fijiq 2 ) and 
one spin-independent Ts{q 2 ) form factors as follows [|38| 



da 
dq 2 



v,q 



8G, 



a 2 ■ T 



00 V 



+ a ai ■ T\q (q 2 
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+ a\.FW) + cl.A^ltf)). (24) 

The last term corresponding to the spin-independent scalar interaction gains 
coherent enhancement A 2 (A is the atomic weight of the nucleus in the reac- 
tion). The coefficients ao,i' c o do not depend on nuclear structure and relate 
to the parameters A q , C q of the effective Lagrangian (pID and to the parame- 
ters Aq, f s , f characterizing the nucleon structure. One has the relationships 

a = {A u + A d ) (An + Ad) + 2AsA s , 

a x = (A u -A d ){Au-Ad), (25) 

Co = f muCu ± ^ + fsCs + hi -fs~ f)(C c + C b + Ct). 
m u + m d 27 

Here Aq p( - n ^ are the fractions of the proton (neutron) spin carried by the quark 
q. The standard definition is 

< p(n)\qYl 5 q\p(n) >= 2S^ n) Aq^ n \ (26) 

where S£, n \ = (0, S p ( n )) is the 4-spin of the nucleon. The parameters Aq p ^> 
can be extracted from data on polarized nucleon structure functions |39], 
and hyperon semileptonic decay data 



We use in the analysis Aq values extracted both from the EMC 39] and 



from SMC g0| data. 

The other nuclear structure parameters f s and / in formula (|25|) are 
defined as follows: 

< p(n) | (m u + m d ) (uu + dd) \p(n) > = 2/M p(n) M, (27) 
< p(n)\m s ss\p(n) > = f s M p[n) ^. 

The values extracted from the data under certain theoretical assumptions 
are |H3l: 



/ = 0.05 and /, = 0.14. (28) 

The strange quark contribution f s is known to be uncertain to about a factor 
of 2. Therefore we take its value in the analysis within the interval 0.07 < 



f s < o.3 m, n 



The nuclear structure comes into play via the form factors J-"ij(q ' ), ^Fs(q ) 
in Eq. (^4|) . The spin-independent form factor Ts{q 2 ) can be represented as 
the normalized Fourier transform of a spherical nuclear ground state density 
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distribution p(r). In the analysis we use the standard Woods-Saxon inspired 
distribution |44|. It leads to the form factor 

T s {q 2 ) = [ d 3 rp(r)e^ = 3^^ e -^) 2 , (29) 
J qK Q 

where i?o = (R 2 — 5s 2 ) 1//2 and s ~ 1 frn are the radius and the thickness of 
a spherical nuclear surface respectively, j\ is the spherical Bessel function of 
index 1. 

Spin-dependent form factors Tij{q 2 ) are much more nuclear model depen- 
dent quantities. The last few years have seen a noticeable progress in detailed 
nuclear model calculations of these form factors. For many nuclei of interest 
in DM search they have been calculated within the conventional shell model 



4q| and within an approach based on the theory of finite Fermi systems ||6 
We use the simple parameterization of the q 2 dependence of Tij(q 2 ) in the 
form of a Gaussian with the r.m.s. spin radius of the nucleus calculated in 
the harmonic well potential 1 . For our purposes this semi-empirical scheme 
is sufficient. 

An experimentally observable quantity is the differential event rate per 
unit mass of the target material 



dR 

IK 





I 


1 m x~ 


Jv 



dvf(v)v^(v,E r ) (30) 



Here f(v) is the velocity distribution of neutralinos in the earth's frame 
which is usually assumed to be a Maxwellian distribution in the galactic 
frame. v max = v esc ~ 600 km/s and p x = 0.3 GeV-cm~ 3 are the escape 
velocity and the mass density of the relic neutralinos in the solar vicinity; 
Vmin = (MAE r /2M 2 ed ) 1 ^ 2 with Ma and M re d being the mass of nucleus A 
and the reduced mass of the neutralino-nucleus system, respectively. Note 
that q 2 = 2M A E r . 

The differential event rate is the most appropriate quantity for compar- 
ing with the observed recoil spectrum and allows one to take properly into 
account spectral characteristics of a specific detector and to separate the 
background. However, in many cases the total event rate R integrated over 
the whole kinematical domain of the recoil energy is sufficient. It is widely 
employed in theoretical papers for estimating the prospects for DM detection, 
ignoring experimental complications which may occur on the way. Notice, 
that the integrated event rate is less sensitive to details of nuclear structure 
then the differential one (|30D . The q 2 shape of the form factors Fijiq 2 ), Fs(<l 2 ) 
in Eq. (^4|) may essentially change from one nuclear model to another. In- 
tegration over q 2 as in the case of the total event rate R reduces this model 
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dependence. In the present paper we are going to perform a general analysis 
aimed at searching for domains with extraordinary large values of the event 
rate R like those reported in This is the reason why we use in the 

analysis the total event rate R. 



5 Numerical Analysis 

In our numerical analysis we randomly scan the MSSM parameter space 
within a broad domain 

1 GeV < mi/2 < 5 TeV, < 2 TeV, (31) 

1 < tan/5 < 50, \Ao\ < 1 TeV, (32) 

< m < 5 TeV, 50 GeV < m A < 1 TeV (33) 

In the region where tan/3 > 35 the top Yukawa dominance approximation 
is not applicable in the RGE. Therefore, we use the procedure developed in 
[O] which takes into account the bottom and tau Yukawa couplings as well. 

The cut-off condition R > 0.01 event/kg/day is implemented in the scan- 
ning procedure. It reflects realistic sensitivities of the present and the near- 
future DM detectors. 

Note again, that we use the GUT scenario with the non-universal Higgs 
mass parameters (see ([J)) Wtfl - The Higgs boson masses are calculated in 



terms of the CP-odd Higgs boson mass rriA and other input parameters. If 
we adopt the ultimate GUT conditions with all scalar mass parameters being 
equal at the unification scale Mx, the CP-even Higgs boson mass becomes 
too big because of strong correlations with the sfermion spectrum. As a 
result the total event rate R decreases to small values, typically less than 
0.01 event /kg/day. 

The main results of our scan are presented in Figs. 1-4 in the form of 
scatter plots. Given in Figs. 1-3 are the total event rates R for 73 Ge, Al 2 03, 
and Nal versus neutralino mass m x , as well as R versus the ratios R s d/R 
of the corresponding spin- dependent (R s d) part of R to R (R = R s d + R S i)- 
Fig. 4 presents the neutralino relic density Sl x h? as a function of m x . All 
quantities are given with and without the b — > s'j constraint. 

We find that the b — > s'j constraint strongly reduces the MSSM parameter 
space. The restriction leaves about 25% of the points of the MSSM parameter 
space which successfully have passed all other constraints. 

This constraint disfavors negative values (in our notation) of the Higgs 
mixing parameter /i. It also shrinks the allowed domain for the parameter 
77ii/2 an d consequently reduces the allowed domain for the LSP mass m x . 
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As seen from Figs. 1-3 the b — > 57 constraint strongly suppresses those 
points in the parameter space which correspond to the spin-dominant (R s d > 
R S i) event rates in all isotopes analyzed. The observation strengthens the 
conclusion about dominance of the spin-independent neutralino interaction 



with nuclei obtained in jLl]] without b — ► S7 constraint. 

Nevertheless it is clear that the large event rates survive the b — > S7 
constraint. In the table 1 we present 5 examples of large event rate points 
taken from the scatter plots in Fig. 1. 

In paper ]10| extraordinary large event rates of about 10 events/kg/day for 
76 Ge and 100 events/kg/day for Nal were found in a specially arranged scan in 
the domain 800 GeV < m x < 1200 GeV, 0.01< Z g <0.99, 0< m A <60 GeV 
(Z g = + M11 see Eq. fll2D). We have thoroughly scanned this region to 
check the cited striking result. We arrived at a negative conclusion. No large 



event rate domains around m x ~ 1 TeV as quoted in |I| have been found 
in our scan. Note, since the neutralino is the LSP, these domains correspond 
to a situation when all SUSY particles are very heavy with masses around 
1 TeV or larger. Looking at the formulas (|22| ) and (|23| ) we do not see any 
natural possibility for R to approach such large values in this domain. The 



strong kinematical suppression can only be compensated in the case when 
rriq — m x ~ m q . 



6 Conclusion 

We have systematically studied the allowed MSSM parameter space taking 
into account various theoretical and experimental constraints. We have found 
domains with experimentally interesting event rates for the DM neutralino 
detection (R ~ 10 events/kg/day) in the neutralino mass range 70 GeV 
< m x < 200 GeV. This would be within the reach of current dark matter 
experiments. Special attention was paid to the constraint following from the 
CLEO measurement of BR(b —>■ 37) . We have illustrated that despite the 
well known fact that this constraint essentially reduces the allowed MSSM 
parameter space it does not exclude large event rate domains. We have 



checked the recently reported result |10| on large neutralino detection event 



rates in the 1 TeV region of the neutralino mass. Our analysis has not 
reproduced this result. 
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Figure legends 



Figure 1. The total event rate R for 73 Ge, versus mass of neutralino m x (upper 
panel) as well as versus the ratio R s d/R of the corresponding spin- 
dependent R s d part of R to the R (R = R s d + R S i)- The scatter plots 
are obtained without (left panel) and with the b — > 57 constraint (right 
panel). 

Figure 2. The same as in Figure 1, but for sapphire, Al 2 03. 

Figure 3. The same as in Figure 1, but for sodium iodide, Nal. 

Figure 4. The relic density £l x h 2 versus mass of neutralino m x without (left panel) 
and with the b — > constraint (right panel). 
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Table 



SUSY points 


1 


2 


3 


4 


5 


tan p 


20.4 


21.2 


21.2 


12.7 


19.5 


m (GeV) 


3654 


1421 


3055 


646 


590 


m 1 / 2 (GeV) 


621 


229 


405 


372 


320 


A (GeV) 


-2.8 


-0.18 


0.5 


-5.3 


-1.4 


m A (GeV) 


941 


588 


673 


575 


685 


H (GeV) 


176 


575 


678 


606 


652 


m x (GeV) 


70.1 


91.3 


163 


148 


129 



0.17 0.074 0.1 0.054 0.16 



R(events/kg/ day) 



7.38 1.08 2.13 2.25 1.73 



R s dl R S i • io 3 



0.5 



0.2 0.04 0.07 0.06 



Gaugino fraction, Z g 



0.04 0.95 0.95 0.93 0.96 



BR(b -> try) TO 3 



0.28 0.22 0.29 0.25 0.19 



Table 1. Representative points with large event rate values for germanium, 
73 Ge. (The gaugino fraction is defined as Z g = Af^ + Af? 2 ) 
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